First-principles investigation of Ag-Cu alloy surfaces in an oxidizing environment 
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In this paper we investigate by means of first-principles density functional theory calculations the 
(111) surface of the Ag-Cu alloy under varying conditions of pressure of the surrounding oxygen 
atmosphere and temperature. This alloy has been recently proposed as a catalyst with improved 
selectivity for ethylene epoxidation with respect to pure silver, the catalyst commonly used in 
industrial applications. Here we show that the presence of oxygen leads to copper segregation 
to the surface. Considering the surface free energy as a function of the surface composition, we 
construct the convex hull to investigate the stability of various surface structures. By including the 
dependence of the free surface energy on the oxygen chemical potential, we are able compute the 
phase diagram of the alloy as a function of temperature, pressure and surface composition. We find 
that, at temperature and pressure typically used in ethylene epoxidation, a number of structures 
can be present on the surface of the alloy, including clean Ag(lll), thin layers of copper oxide and 
thick oxide-like structures. These results are consistent with, and help explain, recent experimental 
results. 

PACS numbers: 68.43.Bc, 81.65.Mq, 61.82.Bg 



T3 

C 

O 

o 



> 
in 
m 

o 



X 



I. INTRODUCTION 



Silver is the active element in heterogeneous cataly- 
sis for a number of industrially important chemical re- 
actions including ethylene epoxidation (carried out at 
atmospheric pressure and temperatures T=500-600 K) 1 
and partial oxidation of methanol to formaldehyde (con- 
ducted at atmospheric pressure and T=800-900 K). The 
origin of the high catalytic activity of such a noble 
metal has been attributed to the thin oxide-like struc- 
tures present on the lowest energy surface for realistic 
operating conditionsi^i^ In agreement with the Sabatier 
principle (in which a good catalyst readily dissociates ad- 
particles but does not bind the fragments too strongly), 
it has been argued that while clean Ag(lll) binds the ad- 
sorbates too weakly, the thin oxides provide O-Ag bonds 
of intermediate strength. A number of energetically very 
similar oxide-like structures have been identified for this 
system 3 -, suggesting a scenario in which the surface might 
comprise of not simply a single low energy structure, but 
could well be dynamically evolving in time, fluctuating 
between different structures. This picture also emerges 
from recent studies on CO oxidation at RuO2(110)£ and 
Pd(100) 6 : In both cases the highest catalytic activity 
is reached in regions of the phase diagram correspond- 
ing to boundaries between different stable surface oxide 
structures. Experiments conducted under semi-realistic 
catalytic conditions for CO oxidation at Pt(110) and 
Pd(100}£ have also shown how, under steady-state catal- 
ysis, the thin oxides present at the surface continuously 
evolve with time. These recent results have therefore 
stressed how in catalysis, a non equilibrium process, the 
catalyst surface cannot be viewed as a static object but 
rather as a "living" system that, at the atomistic level, 



continuously evolves due to various processes such as ad- 
sorption, desorption, association, dissociation and diffu- 
sion. 

Having a detailed knowledge of the catalyst surface 
structures under operating conditions is therefore a cru- 
cial first step toward understanding the full catalytic 
process at the atomistic level. Conventional surface sci- 
ence investigations are conducted in ultrahigh vacuum 
(UHV) and at room temperature or below, while indus- 
trial applications of heterogeneous catalysis usually re- 
quire pressure of the order of atmospheres and tempera- 
tures often higher than room temperature. The informa- 
tion extracted from such surface science experiments can 
not always be extrapolated to realistic conditions, since 
structures that can exist at high temperature and pres- 
sure might not be seen in UHV conditions and viceversa. 
Bridging this difference in ambient conditions (usually 
referred to as the pressure and temperature gap) is one 
of the major goals in current surface science research. To 
this end, first-principles simulations that combine accu- 
rate electronic structure methods with equilibrium ther- 
modynamics have been successfully employed to predict 
the stability of surface structures of several systems under 
high temperature and high pressure conditions^ More- 
over, such methods have been included in a multiscale 
approach that employs statistical mechanics methods to 
simulate the full steady-state catalysis cycle. In some 
cases the agreement with experiments can be quantita- 
tive, allowing an unprecedented level of insight into the 
concerted actions of a dynamical process such as hetero- 
geneous catalysis £ 

With regard to the silver catalyst, the mechanism 
of ethylene epoxidation has been studied in two recent 
works 9 ''' using a combination of experimental and the- 
oretical techniques. It has been proposed that both the 
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selective (producing ethylene oxide) and unselective (pro- 
ducing the acetaldehyde intermediate and finally leading 
to total oxidation) reaction pathways have a common in- 
termediate, a surface oxametallacycle. This intermedi- 
ate can react to form either ethylene oxide or acetalde- 
hyde with similar activation barriers. The same authors 
also reported, on the basis of both first-principles calcula- 
tions^ and experiments^ 2 - that if an Ag-Cu alloy, rather 
than pure Ag, is used as a catalyst, the selectivity toward 
ethylene oxide will be improved. 

The use of bimetallic catalysts like Ag-Cu has been the 
focus of much work in the field of heterogeneous cataly- 
sis^ since the catalytic activity and selectivity of a metal 
can be modified substantially by alloying with another 
metal. For example, geometrical and electronic effects 
obtained by varying the alloy composition may play a 
crucial role in determining the properties of the cata- 
lyst.— One appealing aspect of bimetallic alloys is the 
possibility of rationally designing the catalytic properties 
of the material by changing the alloy elements and com- 
position. To this end, first-principles calculations have 
been shown to be a potentially useful tool for screening 
among pools of possible alloys and for extracting trends 
and therefore gaining insights into the functioning of the 
alloy 15 i 16 . In view of the above discussion, though, the 
effects of high pressure and temperature and the role of 
the dynamical evolution of the catalyst must be carefully 
accounted for, in order for the theoretical simulations to 
have predictive power. 

Through ex-situ X-ray photoelectron spectroscopy 
(XPS) measurements of Ag-Cu catalysts, it has been 
shown that the copper surface content is much higher 
than the overall copper content of the alloyjii therefore 
suggesting copper segregation to the surface. This lead 
Linic et a/£ilil to theoretically model the surface of the 
Ag-Cu alloy assuming a perfect Ag(lll) surface in which 
one out of four silver atoms is replaced by a copper atom. 
However, it is known that copper, at the temperature T 
and pressure po 2 used in such experiments (T=528 K, 
Po 2 =0-1 atm), oxidizes to CuO^i while at higher tem- 
perature or lower pressure CU2O is the stable oxide^ 
Thus it is possible that more complex structures involv- 
ing copper oxide can be present on the catalyst surface. 

The aim of this paper is to investigate, by means of 
first-principles electronic structure calculations, the sur- 
face structure of such alloys under varying oxygen pres- 
sure and temperature conditions. To do this we consider 
the alloy surface to be in thermodynamic equilibrium 
with an atmosphere of pure oxygen. The effect of the 
presence of other reactants such as ethylene is not in- 
vestigated in this work. Under the conditions of steady- 
state catalysis it is therefore likely that the surface of 
the catalyst might be modified with respect to the stable 
structures found in an oxygen atmosphere. Hence our 
work can be thought of as a first step in gaining some 
insight in the possible structures of the alloy, an essential 
prerequisite for modeling the full catalytic process. 

The paper is organized as follows: In Sec.|H]we briefly 



review the theoretical background relevant for this work, 
present the definitions of the quantities that will be used 
throughout the paper and discuss the approximations 
and assumptions that underpin our methodology. In 
Sec. IIIII we report our results and discuss their relevance. 
Section [IVI summarizes the main findings. 

II. CALCULATION METHOD 

The density functional theory (DFT) calculations pre- 
sented in this work are performed using the general- 
ized gradient approximation (GGA) of Perdew-Burke- 
Ernzerhof (PBE)i£ for the exchange and correlation func- 
tional. We use ultrasoft pseudopotential a 19 ' 20 for the 
electron-ion interactions, including scalar relativistic ef- 
fects. The Kohn-Sham wave functions are expanded in 
plane waves with an energy cutoff of 27 Ry (200 Ry 
for the charge density cutoff). To sample the Brilloum- 
zone we use the special-point technique,-- broadening the 
Fermi surface according to the Marzari-Vanderbilt cold- 
smearing technique,- 2 using a smearing parameter of 0.03 
Ry. In the (lxl) surface unit cell, corresponding to the 
periodicity of clean Ag(lll), a 12x12x1 k-point mesh is 
used: This amounts to 19 k-points in the irreducible part 
of the Brillouin-zone and gives adsorption energies for 
oxygen adsorption converged to within 16 meV— at full 
coverage. For larger unit cells, the k-point mesh is scaled 
accordingly. All the calculations are performed using 
the PWscf code contained in the Quantum-ESPRESSO 
packaged For the calculations of chemisorbed structures 
we employ a 7 layer symmetric slab geometry, with ad- 
sorbates on both sides of the slab and fixing the position 
of the three central layers. To model oxide-like structures 
we place the thin copper layers on one side of a three layer 
Ag slab, fixing the bottom two layers 2 ^. The positions 
of the other atoms are relaxed until the forces are less 
than 0.001 Ry/au (0.025 eV/A). A 12 A vacuum layer 
is used, which is found to be sufficient to ensure negligi- 
ble coupling between periodic replicas of the slabi^ For 
the in-plane lattice spacing we use the calculated equilib- 
rium bulk fee lattice parameter of Ag, ao = 4.16 A. The 
experimental value for the lattice parameter is 4.09 Ai^i 

We now define some quantities that will help the dis- 
cussion presented in the following sections. In the case of 
adsorption of oxygen on a clean Ag surface, the average 
binding energy per oxygen atom is defined as 

E° /Ag = L[#0/Ag _ (£; slab + NqE 0^ f (1) 

Ao 

where Ao is the number of oxygen atoms in the unit cell, 
are the energies of the total sys- 
tem, the Ag slab and half the oxygen molecule (which 
has been computed through a spin-polarized DFT calcu- 
lation), i.e. E° = l/2i?Q 2 tal . Defined this way, a positive 
binding energy means that the adsorption of an oxygen 
atom is exothermic (i.e. stable) with respect to oxygen 
in molecular gas phase form. We point out here that the 
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significant error (~ 0.54 eV) in the oxygen molecule bind- 
ing energy introduced by the GGA-PBE approximation 
(1/2E® 2 — 3.10 eV, while the experimental value is 2.56 
eV— ) would introduce a large error bar in the determi- 
nation of the oxygen atom energy and therefore in the 
oxygen atom chemical potential. This error, however, 
will partially be compensated by the analogous GGA- 
PBE error in the description of oxygen chemisorbed on 
the metal surfaces. 

To take into account the effects of temperature (T) 
and pressure (p) we employ "ab initio atomistic ther- 
modynamics" ) 29 i 30 i 31 which allows the determination of 
the lowest energy structures as a function of the T and p. 
The surface is considered to be in contact with an oxygen 
atmosphere that acts as a reservoir, therefore exchang- 
ing oxygen atoms with the surface without changing its 
temperature and pressure (i.e. its chemical potential). 

The change in Gibbs free energy is calculated as 
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where G°' Ag is the free energy of the adsor- 
bate/substrate structure, G slab is the free energy of the 
clean Ag slab, AA^Ag is the difference in the number of 
Ag atoms between the adsorption system and the clean 
Ag slab, /XA g and no are the atom chemical potentials of 
Ag and O. The change in Gibbs free energy is normal- 
ized by the surface area A to allow comparisons between 
structures with different unit cells; we will refer to this 
quantity as "change in Gibbs surface free energy of ad- 
sorption" , or simply "surface free energy" . The chemical 
potential of Ag is taken to be that of an Ag atom in 
bulk Ag, therefore assuming that the slab is in equilib- 
rium with the bulk, that acts as the silver reservoir. The 
oxygen chemical potential depends on temperature and 
pressure according to 3 -^ 
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Ho(T,p) = -[E^ a \T,p )+no 2 (T,P )+kBTln 

2 2 \ p' 

(3) 

Here p° is the standard pressure and /io 2 (T,p°) is the 
chemical potential at the standard pressure, which can 
obtained either from thermochemical tables 3 - 3 , (the choice 
made in this study) or directly computed. 

The vibrational contributions to the free energy, which 
in principle should be accounted for, have been shown to 
be sufficiently small so as not to play an important role 
for the O / Ag system i 30 i 31 When comparing systems with 
different stoichiometry, on the other hand, vibrational 
contributions might play a non-negligible role. In this 
work however we will neglect such effects. To show that 
this approximation is valid, we have follow the procedure 
described in Ref^ to estimate the vibrational contribu- 
tion to the surface free energy, using the Einstein model 
and approximating the phonon density of states by just 
one characteristic frequency for oxygen. We find that for 
two of the most stable structures ("p2" and "CuO(lL)", 
see Sec. IIII[) the oxygen characteristic frequencies are 62 



and 69 meV respectively, giving rise to differences in free 
surface energy of less than 10 meV/A 2 in the range of 
temperatures of interest (0-1000 K). Hence, within this 
model, the only term that depends on T and p is the 
oxygen chemical potential. As a consequence, the free 
energies G°/ Ag and G slab are identified as the total en- 
ergies £°/ Ag and E slah . 

If we now consider the presence of copper impurities 
in the silver surface, we need to modify the definition of 
surface free energy. As we will show later in the paper, we 
will deal with structures in which an overlayer containing 
O, Ag and Cu is adsorbed on a clean Ag surface. If wc 
consider Cu to be in equilibrium with a bulk Cu reservoir, 
the definition of surface free energy becomes: 



AG( MC 



2_( G 0/Cu/Ag 



G 



slab 



- A N Ag a« A g - A Cu ncu - N Mo ) 
(4) 

where Nq u is the number of Cu atoms and /iCu the copper 
chemical potential. Since in this work, as we will show 
in Sec lIHi we are interested in identifying the structures 
belonging to the convex hull of the free energy vs. copper 
content curve, which depends on the curvature of such 
curve, the choice of the Cu chemical potential is arbitrary, 
since the surface free energy depends linearly on it (see 
Eq. 3]) . We now define the average oxygen binding energy 
in the whole structure as 
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(5) 



Using the above-mentioned approximations, we can ex- 
press the change in Gibbs surface free energy in terms of 
the average oxygen binding energy: 



AG(A M o) = -i(iV ^ 6 0/Cu/As + iVoAAto) ,(6) 



where the oxygen chemical potential is now measured 
with respect to half an isolated O2 molecule: A/iq = 

, , 1 ptotal 

Mo - 2^02 ■ 

In writing the change in free energy as in Eq. ([6]) we 
assume that the configurational entropy contribution to 
the free energy is negligible. To show that this is indeed 
the case, we compute the mixing entropy 3 -^ 



S mi * = -k B [x a s ln{x a s ) + (1 - a£)in(l - x a s )] 



(7) 



where is the concentration of the disordered phase a 
and fee is the Boltzmann constant. The quantity — TS nllx 
around temperatures of interest for catalysis (T=600 K) 
has a maximum for x"=0.5, where it takes a value of 
5 meV/A 2 . As we will show later, such contribution is 
negligible compared to the energy scale of the surface free 
energies of the structures considered in this work. 
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III. RESULTS 

A. Copper segregation to the surface 

Copper and silver are known to be almost completely 
immiscible in the bulk due to their large size mismatch 
(13%): The maximum solubility limit is 8.2 mass % Cu at 
the eutectic temperature of 779 °C, and rapidly decreases 
at lower temperatures, with a value of 0.7 mass % Cu at 
400 °C— Accordingly, a DFT calculation at K shows 
that the alloy mixing energy, defined as 

^mix = (£ ,total _ NAgfiAs _ A CuMCu )/(A Ag + A Cu ), (8) 

where E total is the total energy of the unit cell of the 
alloy structure, is positive for all the concentrations con- 
sidered here (between 1.4 % and 6.2 %). The results are 
shown in Table HI Extrapolating the results at zero Cu 
concentration gives a positive mixing energy, therefore 
inhibiting alloying at K at all concentrations. Here we 
relax all the atomic positions and keep the simulation cell 
fixed and commensurate to the computed equilibrium Ag 
lattice constant. One would therefore expect Ag and Cu 
to phase segregate at low temperatures, where entropic 
effects are negligible. 



cell 


%Cu 


iVcu 




cr m,x (meV/atom) 


2x2x1 


6.2 


1 


16 


18.4 


2x2x3 


4.2 


2 


48 


12.6 


2x2x2 


3.1 


1 


32 


9.5 


2x2x3 


2.1 


1 


48 


6.5 


2x3x3 


1.4 


1 


72 


4.4 



TABLE I: Alloy mixing energy cr mlx as a function of the con- 
centration of Cu impurities (in percent) in bulk Ag. The 
simulation cells used are indicated in the first column, given 
in terms of the conventional unit cell vectors, which have been 
kept fixed at the equilibrium Ag lattice constant. Atot is the 
total number of atoms and Ac u is the number of copper atoms 
present in the simulation cell. 

If rather than bulk solid solutions we consider the pres- 
ence of surfaces, we have to note that the surface en- 
ergy of the most stable surfaces of Ag and Cu differ 
significantly: The computed surface energy of Ag(lll), 
7 A g (iii) = o 047 eV /A 2 (exp^ 0.078 eV/A 2 ), is much 
smaller than that of Cu(lll), 7 Cu ( m ) = 0.076 eV/A 2 
(exp^l 0.114 eV/A 2 ). 

This suggests that, when copper impurities are intro- 
duced in silver, it is unlikely for copper to be exposed on 
the surface. To verify this, we compute the total energies 
of Ag slabs in which one substitutional copper atom is 
positioned in different layers, and all the atomic coordi- 
nates are fully relaxed, with the in-plane dimensions of 
the simulation cell fixed to the Ag lattice constant. We 
find that both bulk and surface positions are unfavored 
compared to the position directly below the first layer. 
These results are summarized in Table |TTJ We can see 




FIG. 1: (color online) The (3x3) unit cell used for the study 
of segregation properties of Cu impurities in bulk Ag. The 
dark (blue) sphere represents a Cu atom. 

that by increasing the size of the surface cell the quali- 
tative picture is unchanged. In Fig. [T] we show the unit 
cell used in the case of the (3 x 3) cell. 



Cu @ A-th Ag-layer I AS(meV) 
(lxl) (2x2) (3x3) 

1 56 93 76 

2 -200 -68 -74 

3 -103 -38 -39 

4 



TABLE II: Difference in total energy between the reference 
structure (where Cu occupies an Ag site in the center of the 7 
layer slab, i.e. in layer 4) and with the Cu impurity positioned 
in other layers. The Ag-layer in which Cu sits is indicated in 
the first column. The (lxl), (2x2) and (3x3) surface unit 
cells have been considered. 

Accordingly, we also find that increasing the content 
of Cu in the first layer of the slab leads to an increase of 
the surface energy. This is shown in Fig. [2] as the curve 
with black dots. The positive slope of the curve indicates 
that it is unfavorable for Cu impurities to migrate to the 
surface with respect to staying under the surface or in 
the bulk. 38 The fact that the curve is very slightly con- 
vex (i.e. the black dots lie below the straight line joining 
the two ends at composition and 1), on the other hand, 
indicates that alloying in the first layer is slightly favor- 
able with respect to phase separation into two structures 
that have on the first layer pure Ag and pure Cu. We 
note that on the scale of Fig. [5] this is difficult to see, 
but is evident when plotted on a smaller scale. This be- 
havior (slightly convex curve) is what has been found in 
the theoretical work by Christensen et al£& for the Ag- 
Cu surface alloy on the Cu(100) surface, employing the 
linear muffin-tin orbitals (LMTO) method in the tight- 
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FIG. 2: (color online) Change in Gibbs surface free energy at 
T = K (cf. Eq. HJJ as a function of Cu concentration in the 
first layer of Ag(lll). The oxygen coverages considered are 
indicated in the legend in monolayers (ML). The energy zero 
corresponds the surface energy of the clean Ag(lll) surface. 
The straight dashed lines are obtained by joining the points 
corresponding to pure Ag and pure Cu in the first layer. 

binding representation, using the atomic-sphere approxi- 
mation (ASA) and the coherent-potential approximation 
(CPA) and the LDA approximation for the exchange and 
correlation^ 

The small tendency to alloy and the large difference 
in surface energy between Ag(lll) and Cu(lll) have 
dramatic structural effects as has been recently shown 
from ultra high vacuum (UHV) experiments in which 
Cu thin-film growth on Ag(lll) was studiedi 39 i 40 It was 
found that at room temperature, upon deposition of Cu 
on Ag(lll), Cu forms islands that are encapsulated by 
one monolayer of Ag, in agreement with our DFT results 
presented in this work (cf. Table HI]). 

When oxygen is introduced into the picture, the sit- 
uation changes completely: We find that the presence 
of oxygen chemisorbed on the alloy surface has the re- 
markable effect of reversing the slope of the (black dots) 
curve for the surface energy versus Cu surface composi- 
tion, shown in Fig. [21 For these calculations we consider 
an Ag-Cu alloy in the first layer of the Ag slab, with oxy- 
gen adsorbed in the fee site, i.e. in the most favorable 
adsorption site. Out of the four possible fee sites in the 
(2x2) cell, we always select the ones next to the copper 
atoms present on the first layer for the oxygen adsorption. 
Here we modelled the surface with a 7 layer slab with al- 
loy structures created on both sides of the slab. By fully 
relaxing the structure, we calculate the change in surface 
free energy at T = K as a function of the Cu content 
in the first layer and the oxygen coverage. The results 
are shown in Fig. [5] We can see that the presence of at 
least a quarter of a monolayer (ML) of oxygen induces 
Cu to segregate to the surface, i.e. changes the slope of 



the curve from positive (for an oxygen coverage relative 
to the underlying Ag(lll) lattice O < 0.25 ML) to neg- 
ative (for #o ^ 0.25 ML). The driving mechanism here 
is the strong affinity between oxygen and copper, which 
more than compensates the unfavorable surface energy 
of Cu with respect to Ag: The adsorption energy (cf. 
Eq. (JTJ) of oxygen at 0.25 ML coverage on Ag(lll) and 
Cu(lll) is 0.38 eV/atom and 1.57 eV/atom, respectively. 
We also see that in the absence of copper (i.e. the left 
end of the graph in Fig. the most stable structure is 
the one with Oo = 0.25 ML. At higher oxygen coverages 
the oxygen-oxygen repulsion overcomes the formation of 
O-Ag bonds, leading to an increase in surface energy with 
the oxygen coverage. At the opposite end of the graph, 
where a full monolayer of copper is present in the first 
layer of the slab, the surface energy decreases with the 
oxygen coverage: This is due to the contribution to the 
surface energy of the formation of strong O-Cu bonds. 



B. Surface alloy and surface oxide-like structures 

Having seen that oxygen induces copper to segregate 
to the surface, we now consider various surface structures 
with different contents of copper and oxygen in the first 
layer. These structures, formed on top of a pure silver 
slab, are periodic, and the largest surface unit cell is the 
(4x4). As we will see in the next section, by changing 
the oxygen chemical potential according to Eq. (|3J) we 
investigate the stability of these surfaces by evaluating 
the surface free energy at various temperatures and pres- 
sures. 

We consider three types of structures: (i) chemisorbed 
oxygen on the Ag-Cu alloy in the first layer of the Ag 
slab (as described in the previous section), (ii) structures 
derived from copper (I) oxide CU2O, whose structure can 
be visualized as trilayers of O-Cu-0 piled up on top of 
each other and (Hi) structures derived from copper(II) 
oxide CuO. 

The first set of structures are the ones considered in 
Fig. [21 and we label them OxML/CuyML, where x and 
y are the content of O and Cu, expressed in monolayers 
with respect to the underlying Ag(lll) surface. Some 
of these structures, in particular the ones with high Cu 
content, involve large atomic rearrangements that are due 
to the large Cu-Cu distances induced by the underlying 
Ag(lll) lattice. As an example, we consider two struc- 
tures, one at low oxygen coverage (O0.25ML/CulML) 
and one at high oxygen coverage (OlML/CulML). In 
the first case we find that oxygen adsorbs at a very re- 
duced height in the fee hollow site of the (111) surface 
(0.43 A compared to 1.22 A for the 0/Cu(lll) system 
at the same coverage). In the second case, half the 
oxygen atoms penetrate below the copper layer, lead- 
ing to a trilayer-like structure with a (2x1) periodicity. 
In this structure half the oxygen atoms are positioned 
0.95 A above the average plane of Cu atoms and half 
0.86 A below it. The atomic geometry is illustrated in 



FIG. 3: (color online) Top view of four surface structures considered: (a) "p2", (b) "p4-OCu3", (c) "CuO(lL)" and 
(d) "OlML/CulML" (see text). The grey spheres represent the underlying Ag(lll) substrate. Copper atoms are shown 
as large dark (blue) circles, and oxygen atoms are the small dark (red) circles. The black lines represent the surface unit cells. 



Fig.Etd). 

For the second set of Cu2 0-like structures, we use the 
label "p2" or "p4" depending on whether the periodic- 
ity of the structure is (2x2) or (4x4) with respect to 
the clean Ag surface. The atomic geometry of the "p2" 
structure is shown in Fig. [3^a) . Here each Cu is linearly 
bonded with two O and each O is bonded to three Cu, 
in a ring-like structure similar to the ones proposed for 
thin oxide-like layers on Ag(lll) 3 and Cu(lll)£k The 
"p4-OCu3" is a p4 structure in which an OCU3 unit has 
been removed (see Fig.|3[b)). In "p4+Cu" and "p4+0" 
a Cu atom and an O atom, respectively, has been added 
in an fee adsorption site. "Cue" is analogous to the struc- 
ture derived from the recently proposed "Agg" structure 
for the (4x4)-0/Ag(lll) system42^ We note that con- 
straining an O-Cu-0 trilayer of Cu20(lll) to be com- 
mensurate with the p2 (or equivalently a pA) periodicity 
of Ag(lll) results in a compression (7%) of the trilayer 
with respect to its equilibrium lattice constant in bulk 
copper oxide. We also consider thicker films of CU2O (up 
to 5 atomic layers) in order to extrapolate the behavior 



of bulk CU2O. As we will show in the next section, the 
two most relevant thin oxide-like structures are "p2" and 
"j?4-OCu3". For these Cu20-like configurations, the av- 
erage binding energy per oxygen atom is 1.41 eV and 1.37 
eV, respectively, and the Cu-0 bond length is between 
1.84 and 1.85 A in both cases. For comparison, in bulk 
CU2O the computed formation energy per oxygen atom 
is 1.26 eV (experimental value 1.75 eV— ) and the Cu-0 
bond length is 1.88 A. 

For the third set of structures, we consider thin lay- 
ers of CuO-like structures in a (2x2) cell. Forcing the 
oxide layer to match the (2x2) lattice of the underlying 
Ag(lll) surface leads to a compression of the Cu-0 bond 
length of about 3% with respect to the bulk CuO value. 
For the one layer CuO structure (labelled "CuO(lL)", 
shown in Fig.[3]Jc)) we find that the Cu-0 bond length in 
the square planar pattern is 1.91 A and the binding en- 
ergy per oxygen atom is 1.16 eV. We consider a number 
of possible defective CuO-like thin layer structures (not 
shown here), but we find them not to be stable. We then 
consider thicker CuO-like films (up to 5 atomic layers) in 
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FIG. 4: (color online) Change in Gibbs surface free energy as a 
function of the change in oxygen chemical potential A/io • The 
chemical potential of oxygen is measured with respect to half 
the total energy of the free O2 molecule. The vertical dashed 
lines separate the regions of stability of the clean Ag(lll) 
surface, the surface oxides and bulk oxides CU2O and CuO. 

order to extrapolate the behavior of bulk CuO. In this 
case we must bear in mind that bulk CuO is poorly de- 
scribed with DFT-PBE: CuO is a strongly correlated an- 
tiferromagnetic semiconductor, with a monoclinic struc- 
ture (a=4.65 A, 6=3.41 A, c=5.11 A, /3=99.5°), where 
Cu is linearly bonded to 2 O atoms and O is bound to 2 
Cu4£ DFT-PBE, on the other hand, predicts CuO to be a 
metal with an almost orthorthorombic structure (a=4.34 
A, 6=4.01 A, c=5.22 A, /3=92.2°), in which each O is 
tetrahedrally coordinated to 4 Cu and each Cu is bonded 
in a square planar geometry to 4 O. The predictions for 
thick films of CuO must therefore be regarded as qualita- 
tive. The computed formation energy per oxygen atom 
in bulk CuO is 1.23 eV (experimental value 1.63 eV 44 ) 
and the Cu-0 bond length is about 1.97 A (experimental 
value 1.67 A). Going beyond the DFT-PBE description 
is a major task. Calculations for CuO and CU2O using 
GGA+U plus GoWcr^ as well as other methodology 4 -^ 
are in progress. 

In addition to the structures considered in this work, 
it is likely that other structures involving thin oxide-like 
layers with similar surface energies exist for this system. 
Recent works on 0/Ag(lll)M2>i2 and 0/Pd(lll)^ have 
shown that for these systems a multitude of structures 
with similar energetics can be found. 

C. Thermodynamic diagram of the 0/Cu/Ag(lll) 
system 

In Fig. @] we report the change in Gibbs surface free 
energy of the structures considered as a function of the 
change in oxygen chemical potential. We also show at the 



bottom of the plot, for three values of temperature (300, 
600 and 900 K), the value of pressure corresponding to 
the chemical potential shown in the abscissa of the plot. 

As we can see from Eq. ©, the slope of the lines in 
Fig. 2] is proportional to the oxygen coverage, i.e. the 
higher the oxygen content the steeper the line. The 
vertical dashed line at A/io = —1-26 eV represents 
the change in chemical potential above which Cu oxi- 
dizes to CU2O, which, as one can show, corresponds to 
the computed heat of formation of CU2O. The vertical 
dashed line at A/io = —1-23 eV, on the other hand, 
represents the change in chemical potential above which 
Cu oxidizes to CuO. We can therefore see that in the 
range of values of the oxygen chemical potential con- 
sidered (—1.6 < A/io < —0.4 eV) the thermodynami- 
cally stable structures are: Pure Ag for A/io < —1-43 
eV, "p2" and "p4-OCu3" (almost degenerate) between 
-1.43 < A/i < -1.26 eV, bulk copper (I) oxide Cu 2 
for —1.26 < A/io < —1.23 eV and bulk copper(II) oxide 
CuO for A/io > —1-23 eV. It is interesting to note that 
none of the alloyed chemisorbed structures (the struc- 
tures labelled as OirML/CuyML) are thermodynamically 
stable in this range of chemical potential, while there 
is a small region in which two-dimensional surface ox- 
ides are stable. This situation is similar to what has 
been found for the 0/Cu(lll) system^ where it was ar- 
gued that copper oxidation does not proceed via ordered 
chemisorbed structures, at variance with other transition 
metal structures. We also find that none of the O/Ag 
structures, containing no Cu (not shown in Fig. [4]) are 
thermodynamically stable in the range considered here. 
This is not unexpected, given the weaker strength of the 
Ag-0 bond with respect to the Cu-0 one. 



D. Modeling systems of known Cu surface content 

Having identified, among those considered here, the 
thermodynamically most stable ordered structures on an 
infinite Ag(lll) surface, we now consider the situation in 
which we have a well defined Cu surface content. For a 
finite system (e.g. an Ag nanoparticle) of known dimen- 
sions and well defined Cu content, if we assume that all 
the Cu present in a spherical nanoparticle segregates to 
the surface in a site of the Ag fee lattice, given the size of 
nanoparticle and the total Cu content, we can estimate 
the Cu surface content. For example in a 50 nm spher- 
ical nanoparticle with a 1 % Cu content (typical values 
of nanoparticle size and copper content of Ag-Cu alloys 
used in experiments in which the effect of Cu impurities 
was investigated 1 ^) the (average) surface Cu content is 
0.35 ML, where 1 ML corresponds to the first layer of 
the nanoparticle being pure copper. Experiments carried 
out on Ag-Cu alloys with 0.1-1% Cu content at ~ 500 
K have shown that in an oxidizing environment the Cu 
surface content is in the range 0.1-0.75 MLJ£ 

Given the fact that Ag(lll) is the lowest energy sur- 
face for silver^ the typical understanding is that a large 
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FIG. 5: (color online) Construction of the convex hull for the 
case of A/io = —0.80 eV. The energy zero corresponds to the 
clean Ag(lll) surface. 



FIG. 6: Surface phase diagram showing the structures belong- 
ing to the convex hull as a function of the Cu surface content 
and the change in oxygen chemical potential. 



portion of the nanoparticle surface will consist of (111) 
planes. This is not necessarily true at high temperature, 
where the particle can be quite round, especially in the 
case of noble metals, and in the presence of oxide-like 
structures on the surface, that can alter the relative sur- 
face energies compared to the case of pure Ag at zero 
temperature. We will however restrict ourselves to this 
particular surface, and use the results shown in Fig. [4] 
to predict which structures will be present on the sur- 
face as a function of the copper content and the oxygen 
chemical potential. We also exploit the vast literature 
available for the O/Ag system^^ to include the most 
stable structures for the system in the absence of copper. 

As an example, we show explicitly the case for 
A/io = —0.80 eV, corresponding to an oxygen pressure 
of 10~ 3 atm at the temperature of 600 K. Figure [5] shows 
the change in surface free energy at A/io — —0.80 eV of 
all the structures we have considered (a total of 55 struc- 
tures, including the O/Ag ones) as a function of the Cu 
surface content. By constructing the convex hull, i.e. the 
curve obtained by joining those structures that are sta- 
ble with respect to linear combinations of structures at 
other compositions that would yield the same total com- 
position, we identify the structures (indicated by the red 
dots in Fig. [5]) that are stable against phase separation 
into any two other structures. The structures belonging 
to the convex hull, for this particular value of chemical 
potential, are "Ag" (clean Ag(lll) surface), "p4-Cu30", 
"p2" and "CuO(b)" (bulk CuO). Although in Fig. [5] we 
show only the portion of the diagram up to a Cu content 
of 1 ML, we have computed structures including up to 
five CuO layers in order to be able to extrapolate the 
behavior of bulk CuO. 

For a specific Cu content, e.g. 0.25 ML, the convex 
hull plot helps to predict that a mixture of pure Ag and 



"p4-Cu30" (rather than a single ordered structure with 
0.25 ML Cu content) will be present on the Ag(lll) sur- 
face at A/io = —0.80 eV, in a ratio given by the lever 
rule, i.e. by the mass conservation law. By repeating this 
scheme for all the values of the oxygen chemical potential 
in the range considered in this work, we build a phase di- 
agram as a function of the oxygen chemical potential and 
the copper surface content. This is shown in Fig. [6l We 
find that the convex hull, for high Cu content, always 
includes the structure with the largest number of CuO 
layers (or CU2O layers at lower oxygen chemical poten- 
tial), i.e. bulk oxide formation is favored beyond a certain 
Cu surface content. We therefore cut the plot in Fig. [S] 
at a Cu composition of 1.50 ML, where it is understood 
that the rightmost structure is "bulk" oxide. However, 
care must be taken in the meaning given to such "bulk" 
structure in the context of a finite system as the one we 
are modelling here. We interpret the numerical evidence 
of the presence of the bulk structure in the convex hull 
as a tendency for the formation of thick patches of either 
CuO or Cu 2 0. 

In the phase diagram shown in Fig. [5] we have explic- 
itly written the two structures coexisting in a number 
of regions around the values of interest of oxygen chem- 
ical potential and copper content. The labels "Agi. 5 0" 
and "Agi.20 Asym " refer to O/Ag structures identified in 
Ref. [3|- If we focus on the region around the chemical po- 
tential of interest in typical industrial applications (p=l 
atm, T=600 K, corresponding to A/io ~ —0.61 eV) and 
for contents of copper below half a monolayer, we predict 
patches of one layer oxidic structures ("p4-Cu30") to co- 
exist with the clean Ag surface. At higher values of oxy- 
gen chemical potential, on the other hand, O/Ag struc- 
tures can be found in coexistence with the "p4-Cu30" 
structure. For higher copper contents, the "CuO(lL)" 
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and the "p2" structures are present in the phase diagram 
above and below A/iq ~ —0.75 cV, respectively. Finally, 
for even higher copper contents, bulk CuO is predicted 
to form. 

This picture is consistent with recent experiments 
performed on the Ag-Cu system under catalytic condi- 
tions^ In these experiments Ag-Cu nanopowders ( ~ 
100 nm in diameter and 2.5% Cu) are used as the cata- 
lyst for ethylene epoxidation at T — 520 K and p = 0.5 
mbar. Through a combination of in- situ XPS, and Near 
Edge X-ray Absorption Fine Structure (NEXAFS) mea- 
surements, thin layers of CuO are found be to present 
on the surface, while no signs on a surface alloy are de- 
tected. Areas of clean Ag are also present on the surface, 
in agreement with our theoretical calculations. 

These results show that the simple structure adopted 
in Ref. [ll| for the Ag-Cu surface alloy to model theoreti- 
cally the ethylene epoxidation reaction is not stable and 
is significantly different from what we expect to be rele- 
vant under high pressure conditions. The model used in 
Ref. ITT1 assumes a (2x2) periodicity for the Ag(lll) sur- 
face in which one atom out of four is replaced with cop- 
per, and oxygen is chemisorbed on it. Our results suggest 
on the other hand that a model that includes, depend- 
ing on the Cu surface content, clean Ag(lll), patches of 
oxide-like thin layers and thick layers of CuO is a more 
appropriate model for the Ag-Cu catalyst (111) surface. 

IV. CONCLUSIONS 

In summary, through density functional theory calcula- 
tions and employing ab-initio atomistic thermodynamics 



to include the effects of pressure and temperature, we 
have investigated the (111) surface of the Ag-Cu alloy. 
We have shown that, in the absence of oxygen, copper 
impurities in silver prefer to stay directly below the sil- 
ver surface, rather than in the bulk or on the surface. 
The presence of oxygen, on the other hand, has the effect 
to induce copper to segregate to the surface due to the 
strength of the O-Cu bond relative to the O-Ag one. Both 
findings are in agreement with experimental evidence. 
We have investigated structures involving chemisorbed 
oxygen on the alloy surface, as well as oxide-like struc- 
tures on the Ag(lll) surface. Through the construction 
of the oxygen chemical potential dependent convex hull, 
we have been able to identify, as a function of the sur- 
face copper content, the combinations of structures that 
are most stable at that particular temperature and pres- 
sure. In the region of interest, our results suggest that, 
depending on the copper surface content, clean Ag(lll) 
and thin copper oxide-like structures ("p4-Cu30", "p2" 
and "CuO(lL)") can coexist. This model for the (111) 
surface of the Ag-Cu catalyst differs substantially from 
the structures used in earlier works to investigate the ef- 
fect of copper impurities in silver on the mechanism of 
ethylene epoxidation. Our results suggest that to gain 
insight into the full catalytic cycle one should consider 
oxygen species such as those in the proposed structures, 
rather than just oxygen chemisorbed on the the clean 
alloy surface. We envisage that the theoretical method- 
ology reported in the present paper will also be useful for 
characterizing other alloy catalysts under realistic condi- 
tions. 
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